Tunneling dynamics in relativistic and nonrelativistic wave equations 
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We obtain the solution of a relativistic wave equation and compare it with the solution of the 
Schrodinger equation for a source with a sharp onset and excitation frequencies below cut-off. A 
scaling of position and time reduces to a single case all the (below cut-off) nonrelativistic solutions, 
but no such simplification holds for the relativistic equation, so that qualitatively different "shallow" 
and "deep" tunneling regimes may be identified relativistically. The nonrelativistic forerunner at a 
position beyond the penetration length of the asymptotic stationary wave does not tunnel; never- 
theless, it arrives at the traversal (semiclassical or Buttiker-Landauer) time t. The corresponding 
relativistic forerunner is more complex: it oscillates due to the interference between two saddle point 
contributions, and may be characterized by two times for the arrival of the maxima of lower and 
upper envelops. There is in addition an earlier relativistic forerunner, right after the causal front, 
which does tunnel. Within the penetration length, tunneling is more robust for the precursors of 
the relativistic equation. 
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I. INTRODUCTION 

Tunneling is a general feature of wave equations: for 
frequencies below the "cutoff", wavenumbers become 
imaginary and the waves evanescent. Even so, it has 
been mainly discussed for the Schrodinger equation due 
to the shocking contrast between the behavior of classi- 
cal and quantum particles. It has been also tradition- 
ally examined for stationary conditions, until a paper by 
Biittiker and Landauer [| triggered the interest on its 
temporal aspects, again, mostly within the framework of 
the Schrodinger equation. 

Tunneling dynamics may be very counterintuitive 
when compared with the evolution of propagating waves. 
Defining "tunneling times" has been in particular non 
trivial in quantum theory (Ref. is a recent and 
extense multi-author review; for previous reviews see 
H H S H 0) After twenty years of discussions and 
many proposals since []], it is clear that different time 
scales are relevant depending on the experiment, wave 
feature, or observable quantity investigated. There is 
thus not a single tunneling time scale, but several; many 
of them may be obtained, related to each other, and clas- 
sified systematically by means of a theory that exploits 
the non-commutability of the operators implied @. 

One of the most important and commonly found 
scales is the "traversal", "semiclassical" or "Buttiker- 



Landauer" time r. For the dimensionless Schrodinger 
equation for a particle moving in a constant potential, 



1 d 2 ip 



(1) 



(all quantities are dimensionless throughout this work ex- 
cept in the discussion of section ITTT|l . and for a frequency 
ljq < 1, it is given by Q 



x 

k ' 



where 



«o = [2(l-c )] 1 /2 



(2) 



(3) 



in Eq. J2Jl plays the role of a semiclassical "velocity" 
which increases for deeper and deeper tunneling (i.e. for 
smaller ujo). r is the characteristic time for the transition 
from sudden to adiabatic regimes in oscillating barriers, 
and for the spin rotation in a weak magnetic field. A 
recent and surprising finding is that it is also the time of 
arrival of the peak of the forerunner that appears, beyond 
the penetration length of the asymptotic stationary wave, 
x > 1/kq, for a "source" with a sharp onset, 



yj(x = 0,t) = e-^ ot Q(t), 



(4) 



and 
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yb{x > 0,t < 0) = 0. 



(5) 



Q(t) is the Heavisidc function; in the language of electro- 
magnetic wave propagation it is the envelop function of 
the input pulse, whereas wo is the "signal" , "carrier" , or 
"excitation" frequency. More precisely, the arrival time 



2 



is i* = r/3 1 / 2 if the peak of the density \ip(x, t)\ 2 is eval- 
uated with respect to t for x fixed 0, and = r if it is 
evaluated with respect to x for i fixed [Tol | . 

This role of r is indeed surprising because that fore- 
runner is not tunneling beyond 1/kq, where it is predomi- 
nantly composed by over-the-barrier (u> > 1) components 
0. The importance of r is blurred however within the 
penetration shell, x < 1/ko, where the arrival of the tun- 
neling peak obeys a different time scale, inversely pro- 
portional to Kq rather than to Kq 0, 0] . In the small- x 
region the arrival time of the peak t l v as a function of 
x is rather flat in comparison with the (linear in x) ar- 
rival time of the peak in the large-x region. It actually 
forms a basin with a minimum, so that the precursor 
may arrive later at smaller values of x [Til 

El- 

Moreover 

the forerunner does tunnel at small x. This is another 
unexpected result |12j . after so many works where pre- 
cursors had been always related to an above-the-barrier 
passage El El 13 El Surely no attention had been 
paid to the small- a; region because the analytical approx- 
imations based on saddle and pole contributions of the in- 
tegral defining the solution fail, as discussed below. This 
region is the easiest to observe though. 

Our initial motivation to undertake the present work 
was to determine if and how the above recent results 
extend to the reiativistic case. Especifically, our aim is to 
investigate the role played by r, and the characteristics 
of the forerunners at large and small x for under-cut- 
off, unit-step modulated excitations with the reiativistic 
wave equation 



d 2 d 2 
dt 2 dx 2 



i> = 0, 



(0) 



which is much more accessible experimentally (with 
waveguides in evanescent conditions 

El El El) than 

the Schrodinger equation. These are all aspects that pre- 
vious studies of reiativistic Klein-Gordon equations with 
evanescent conditions had not explored [20L l2l| . Some 
preliminary, partial results on the role of the reiativistic 
r may be found in d . 

Note that the stationary equations corresponding to 
Eqs. (JIJ and © are equal and have equal solutions, 
but the two cases are not equivalent in the time domain, 
because the dispersion relations between u> and k are dif- 
ferent, 



k 2 (reiativistic), 
fc 2 /2 (non reiativistic). 



(7) 

(8) 



In fact we have found a wealth of qualitative changes 
with respect to the Schrodinger scenario as the following 
sections will show. One of them is discussed in section 
ITT1 namely, the absence for the reiativistic equation of the 
scaling properties that simplify the nonrelativistic solu- 
tions; this implies that qualitatively different solutions 
are possible for different carrier frequencies in the rei- 
ativistic case. Section IIIII deals with possible physical 
realizations of the formalism. The technical tools used 



are an analytical series expression for the exact solution, 
which is obtained in section lTVI and the asymptotic anal- 
ysis of section|Vl based on previous work by Biittiker and 
Thomas (2^, that will allow to describe and understand 
the wave behaviour at "large-x" in simple terms. For 
small- x this analysis is substituted by numerical explo- 
ration in section IVT1 except in the proximity of the very 
first front or "causal limit" at t = x, where a simple 
analytical approximation is again available. 



II. SCALING PROPERTIES 

The solution of Eqs. @ or Q with the boundary 
condition of Eq. is given in both cases by 



ip(x,t;oj ) 



2vr 



(kt) 



Akx—iut 



uj — cjq + iO ' 



(9) 



as may be checked by substitution. We shall only con- 
sider x > 0. It is easy to see from Eq. © that the wave 
vanishes for t < x. In the reiativistic case k = (lu 2 — l) 1 / 2 
is defined with two branch cuts from to = ±1 downwards 
in the complex cj-plane, as shown in Figure 1, whereas for 
the nonrelativistic equation, k = [2(uj — l)] 1 / 2 is defined 
with only one branch cut downwards from ui = 1. For 
evanescent waves, i.e., when k — in, k > the reiativistic 
wavenumber tends to the nonrelativistic one as lo — > 1, 



(1-c^ 2 ) 1 / 2 



[2(1-^)] 1/2 . 



(10) 



In the nonrelativistic case it is useful to shift the fre- 
quency axis by defining the variable 



CI 



1 



(11) 



(in particular f^o = u>q — 1). Let us also define 
4>(x, t; Slo) = e lt 4>(x,t;ujo). If we introduce a scale fac- 
tor 

fl = aQ', a > 0, (12) 

it follows from Eq. © that 

<j>(x,t;n ) = ())(a 1/2 x,at-no/a), (13) 

which means, in words, that any two solutions with ex- 
citation frequency above cut-off (f^o > 0) are related to 
each other by the scaling law of Eq. i|13fl and similarly, all 
solutions below cut-off (Slo < 0) are also simply related to 
each other by scaling. We can generate all possible solu- 
tions from two of them, one below and one above cut-off. 
Nevertheless, Eq. (|13|l does not hold for the reiativistic 
wave equation because of the different dispersion rela- 
tion, so that solutions for different excitation frequencies 
wo cannot reduce to each other; they may and do change 
qualitatively when sweeping over ujq, as demonstrated in 
Section llVl 
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III. PHYSICAL REALIZATIONS 

Even though we use throughout the paper dimension- 
less variables and equations, it is worth noting that Eq. 
© leads to the Klein-Gordon equation for spin-0 parti- 
cles by means of the substitutions 

x = X(m c)/K, (14) 
t = Tm c 2 /h, (15) 

where X and T denote dimensional position and time, 
and m-o and c are the particle's rest mass and the veloc- 
ity of light in vacuum. Another interesting connection, 
much more promising to implement the present analysis 
experimentally 19] , and free from the conceptual puzzles 
of the former, may be established with the equation that 
governs the electromagnetic field components in waveg- 
uides of general constant cross section with perfectly con- 
ducting walls 23] . This requires the substitutions 

x = XX, (16) 
t = TXc, (17) 




-3-2-10123 



Re(co) 

FIG. 1: Steppest descent paths in the complex w plane for 
t — 1.5, luq — 0.85 and x — 1 (narrow solid line). The thicker 
solid lines are the two branch cuts while the verticals lines are 
the asymptotes. The dots mark ±w s and ±l/w s ; the cross is 
at wo- 



IV. RELATIVISTIC SERIES SOLUTION 



where A > is one of the eigenvalues of the waveguide. 

Dimensional Schrodinger equations corresponding to 
Eq. J3> may be obtained formally in the limit c — > oo 
(the proper one from the Klein-Gordon equation for 
particles, and an analogous one from the waveguide 
equation). Physically, the solution of the dimensionlcss 
wave equation approaches the solution of the dimension- 
less Schrodinger equation when it is dominated by w- 
components close to one (the cut-off), note that the rela- 
tivistic dispersion relation, Eq. J7J), tends to the nonrel- 
ativistic one, Eq. iJHJ, as w — > 1. Indeed excitations with 
wo ~ 1 provide relativistic solutions very close to the non- 
relativistic ones except in the region near the causal limit 
t = x (the "first precursor"), where the whole frequency 
range of the source signal has a significant influence. 

In order to compare relativistic and nonrelativistic 
equations the same "excitation" will be used, Eq. 
which for the waveguide may be understood as a compact 
representation of (real) sine and cosine excitations. Since 
for the relativistic equation the excitations e ±luJf,t <d(t) 
produce at x the complex conjugate wave of each other, 
only the interval < wo < 1 will be examined. Also, 
to carry out the comparison with the Schrodinger equa- 
tion in terms of analogous quantities, we shall study \4i\ 2 
which, in the relativistic-waveguide case may be simply 
regarded as the sum of the squares of the waves result- 
ing from the sine and cosine excitations. Because of its 
meaning in the Schrodinger equation, we shall refer to 
this quantity as the "density" , even in the relativistic 
case. The agreement between the relativistic and non- 
relativistic solutions near the cut-off and far from t = x, 
combined with the scaling property of Eq. (|13fl may serve 
in practice to simulate in a waveguide the solution of the 
Schrodinger equation for any wq in the time domain. 



In this section we shall obtain the solution of Eq. © 
in series form. It will be used to provide exact numerical 
results as well as analytical approximations for x, t > 
and for the sharp onset source of Eq. Q, following refs. 
p4l |25| . We begin by Laplace transforming the equation 
© using the standard definition, 



ip(x;s) = / tp(x,t)e dt. 
Jo 



(18) 



The Laplace transformed solution reads, 

^(x;s) = c ie ~ px , x>0, (19) 

where we have defined p = (s 2 +l) 1 / 2 . The corresponding 
Laplace transform of Eq. (0J, yields, 



iw 



(20) 



By combining Eqs. i|19|) an d l|2U|l evaluated at x = we 
can determine the value of the constant ci, 



ip(x;s) = 



x>0. 



zw 



(21) 



The time dependent solution for x > is obtained by 
performing the inverse Laplace transform of Eq. 121fl . 
using the Bromwich integral formula, 



i>{x, t) 



1 

2tt7 



i/j(x] s)e st ds, t > 0, 



(22) 



7 — too 



where the integration path is taken along a straight line 
s = 7' parallel to the imaginary axis in the complex s- 
plane. The real parameter 7' can be chosen arbitrarily 
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as long as all singularities remain to the left-hand side of 
s = 7'. 

To avoid dealing with the branch points at s = ±« in 
Eq. Q22p. let us introduce the change of variable, — iu = 

~ u)/2, and the integral may be 



(s+p). Thus, p = i(i 



written as 



27+00 



4>{x,t) 



2ni 



(23) 



£7 — OO 



where the new integrand is given by 

(i-« 2 ) 



$0) 



u(u — luo - k )(u — ujo + k ) 
xexp{i[w(a; - t) - U~ X {x + t)}/2}, (24) 



and kg = k(uJo). The integrand has now an essential 
singularity at u = and two simple poles u± = (ujQ±ko). 
For x > t we close the integration path from above, by a 
large semicircle Ti of radius R, forming a closed contour 
C\. The contribution along Ti vanishes as R — > 00, and 
since there are no poles enclosed inside Ci, ip(x,t) = 
for x > t. For the case x < t, we close the integration 
path from below with a large semicircle T2- The closed 
contour C2 contains three small circles Co, C+ and C_ 
enclosing the essential singularity at u = and the simple 
poles at it+ and w_, respectively. Hence, it follows that, 



1 

2?ri 



47 + 00 



27 — 00 



Co C-f 



= 0. (25) 



The integrals corresponding to the contours C+ and C_ 
can be easily evaluated, 



1 



(26) 



The contour integration for Co involves an essential sin- 
gularity at u = 0. Introducing the change of variable 
V = —iut;^ 1 the integral becomes 



1 

2?ri 



$(u)<fu 



_L_ /• 1 (l + ^je"^-"" 1 )/ 2 
2ni J £ 2 u(« — t | +)(t | — V-) 
c' 



where 



v± 
z± 

v 



z±/it 
{u>a ± fc ), 

[(t + .*)/(t-:r)] 1/2 . 



off , 

(27) 

(28) 
(29) 
(30) 
(31) 



The integration is carried out by first separating the in- 
tegrand into partial fractions, and then substituting the 
formula for the Bessel generating function, 



E vnj »(v) + J2(-^) n v- n Uv), (32) 



and the series expansion, 

oc 



(33) 



n=0 



Now we can calculate the resulting integrals by means 
of the residue theorem. For the case of an essential sin- 
gularity, the residue may be determined by computing 
explicitly the coefficient corresponding to w" 1 from the 
series expansion and their products. In that case, equa- 
tion (|2"T|) becomes, 



1 



$(u)du 



-■/ofa)-£(-l) 

n=l 



. Jn(n) 

(v+) n 



n=l 



(34) 



Finally, by feeding the results given by Eqs. (12(i[) and 
(|34|l into Eq. I|25|) . the solution for the internal region 
reads 



tp{x,t) = 



4>{x, ko, t) + ip(x, — ko, i), t>x 



0, t < x, 

with ip(x, ±fe, t) defined as 

*Kx,±ko,t) = eV^o'-^ + ^Mri) 

OO 



(35) 



(36) 



71=0 



It is worthwhile to remark the similarity of the above 
solution with that of Eq. (16) of ref. and Eqs. (15) 
and (16) of ref. From Eq. J22J, 



(37) 



n— n— 1 

Therefore, Eq. JSJJJl may also be written in the form 



1 00 

iKx,±ka,t) = Mv)+J2(-^ n ( iz ^0 n Mv), (38) 



which will be useful for the analysis of the region close 
to the causal limit t = x. 

The above relativistic solutions should be contrasted 
with the solution to the nonrelativistic equation using 
the same initial condition [jj, E| , 

i>(x,t) = ^e-^ 2 /^[w(y k0 ) + w(y^ )} (39) 

where w(y±k ) stands for the u w- function" [2(| with ar- 
gument y±k , 



y±k = ic-^ 4 



1/2 



■ T k t] 



(40) 
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V. SADDLE-POLE APPROXIMATION FOR 
LARGE x 

Following Sommerfeld and Brillouin |23|, who studied 
the propagation of a unit step-function modulated signal 
in a Lorentz medium 1 , the forerunners and other wave 
features may be understood and quantify with asymp- 
totic analysis techniques. Let us deform the integral of 
Eq. © along the steepest descent paths (SDP), see Fig- 
ure ^ defined by 



UJ, 



-(±W R CJ a - 1)(±Wr 



(41) 



(ujr and u>i are the real and imaginary parts of uj), where 
the upper sign is for the positive saddle at uj s and the 
lower sign for the negative saddle 2 at — uj s , 



±w.=± 



t 



(t 2 



2)1/2' 



(42) 



plus a clockwise circle around the pole at ujq that must be 
added to the contour if the pole has been crossed by the 
SDP. The SDP have asymptotes at ujr = ±[uu s ± (uj 2 — 
l) 1 / 2 ] and cross the real axis at ±uj s and iwj 1 , see figure 
n The crossing of the pole occurs when uj~ x = u)q, i.e., 
for fixed x, at time 



x 

Kq 



(1 



2)1/2' 



(43) 



which generalizes the traversal time of Eq. J2J) for the 
relativistic case j2^|. Note that r > x, i.e., the crossing 
occurs always after the arrival of the first front imposed 
by relativistic causality at t = x. One might be tempted 
to believe that the crossing of the pole is associated with 
the arrival of a "monochromatic" wo-front j2^. That 
this is not the case has been pointed out in several works 
on the Schrodinger equation, the reason being the domi- 
nance of over-the-barrier components associated with the 
saddle. The frequency analysis of the relativistic fore- 
runner is more complex as the following discussion will 
demonstrate. 

Applying the standard asymptotic approximation for 
the saddle contributions [3(]], ip takes the form 



lp(x, t) SW lp p + Ipt + 1p 
UJ 



i>7 = * 



2irt 2 ujs T <+>o 



(44) 
(45) 

(46) 



Eq. (|44|l is a good approximation as long as the "width" 
of the saddle is small compared to the "distance" from 



the saddle to the pole |2S 



8uj 2 



:(oj s Tuj ) 2 » 1, 



(47) 



where <p = ujt — kx. This always happens for large enough 
t and fails very close to t — x. Applying this condition 
at t — t we get 

xojqKq » 1 (positive saddle), (48) 
xuj (1 + luq) 2 /kI » 1 (negative saddle), (49) 

to be compared with the simpler nonrelativistic criterion 



9] 



XKq >> 1. 



(50) 



At variance with the nonrelativistic case, the condi- 
tion of Eq. ll4U|l may only worsen with decreasing u , 
whereas the one in Eq. H48|) improves from loq = 1 up to 
the critical value ujq — (1/2) 1 / 2 s» 0.7, but also worsens 
below, ifjf is dominant with respect to as ujq — > 1 
and for large t but the importance of the negative saddle 
increases with decreasing ujq and t, 



t + 7]LUq 



In particular, at t 



\i>7\ 

-- r, 



t 



T]LU 



1 - wg ' 



(51) 



(52) 



To compare the relative importance of the saddle and 
pole contribution we can examine the ratio \ip p /ipf\, 



TZ(x,t) = 



i 2 



-2kqx 



e(t 



t 2 rj 

2~ 



(53) 



TZ(t, x) grows monotonously with t and decreases with 
x. The time scale for the attainment of the station- 
ary regime, or equivalently, the duration of the transient 
regime dominated by the saddles before the pole dom- 
inates, ttr, can be identified as the time when TZ = 1. 
Assuming that t <^tt r , this time is given by 



,2c,i/3 2/3 



(2tt) 1 /3(i_ Wo )2/3 



(54) 



which is an exponentially large quantity as in the 
Schrodinger equation |9j]. Indeed, Eq. (|54|l tends to the 
nonrelativistic result taking u>o — > 1 ■ 



"Densities" 



1 For a modern and more accurate treatment see \>A 

2 Compare with the only nonrelativistic saddle at lu b = 1 + 
x 2 /(2t 2 ), which is obtained from the positive relativistic saddle 
for t >> x. 



If the pole contribution is negligible compared to the 
two saddle point contributions, the total density is given 
by 

hMM)| 2 «|^- + d 2 = 



6 
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FIG. 2: \i>+ 1 2 (solid black line) , \ip~ | 2 (dashed line) and + 
ip7\ 2 (solid grey line); \ips\ 2 (nonrelativistic: dots), x — 100 
and u)o = 0.99. 
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FIG. 3: Upper (a) and lower (b) envelops at at x — 100 for : 
ujo = 0.99 (solid line), ujq — 0.97 (dashed line) and u>o = 0.5 
(long-dashed line). 



t 2 + u%rf + x 2 {t 2 /r 2 + uj 2 ) sin(2r?) 
Trrjx 2 (t 2 /r 2 + uj 2 ) 2 



n < i (55) 



The interference between the two saddles leads to an 
oscillatory pattern of the density, whose amplitude A in- 
creases for deeper tunneling, see Figure [21 From Eq. 



A 



TTTjit 2 



uj 2 r] 2 ) ' 



(56) 



The oscillation only disappears in the nonrelativistic 
limit uj — > 1, t >> x. 

Since the dominant frequencies of the saddle terms for 
tpf and are ±w s respectively, the oscillation period of 



FIG. 4: t^Jr (dashed line); t^Jr (solid line); t^/r 
(dotted-dashed line); tl^ in /r (dotted line); ^^/t (circles), 
and t\£ ax /T (crosses) 



the interference pattern between the two saddles is very 
well approximated by tt/lu s . Of course, beyond t tr the 
oscillations fade away and are substituted by the asymp- 
totic dominance of the pole term. 

are given, 



The upper and lower envelops of \tpf ■ 
from Eq. lJ55j). by: 



2i 2 



*7\l 



■KX 2 r\ (t 2 /t 2 + uj 2 y 
tufa 

TTX 2 (t 2 /T 2 +UJ 2 f 



(57) 



A second forerunner, after the first one at t = x, may thus 
be identified thanks to the maxima of these envelops, see 
Fig. |21 The upper envelop may hold a maximum and 
a minimum, but they disapear below a critical value of 
the injection frequency u>q, see Fig. [3] a. We may thus 
clearly distinguish between shallow and deep tunneling 
regimes with qualitatively different features for the up- 
per envelop. On the other hand, the lower envelop has 
only a maximum that remains for any u>o, although its 
amplitude decreases with smaller ojq, see Eq. (|57|l and 
Fig. 0b, and becomes physically irrelevant in comparison 
to the upper envelop. 

Taking the derivative with respect to t or x of the en- 
velops, we can extract "temporal" or "espacial extrema" 
respectively. In the case of the lower envelop, 



max 



3-u; 2 + ^9(l + <4) 



Uu 2 



1/2 



I ... 



4-3^ 



1/2 



(58) 



which are always larger than the corresponding nonrel- 
ativistic limits, t and r/y/3 respectively. For the upper 
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-0,02 



FIG. 5: \tp\ 2 and u> av versus time for x — 100 and loo = 0.5 
(dashed line and solid line respectivity) . Also shown is the 
nonrelativistic uo a v (dotted-dashed line). 



envelop we obtain 



t*' D , . 

max j min 
,x,U 



max I m%r: 



T 
T 

2 



4 ± y / 4-28wg + 25o;4 



1/2 



1 + ± y 1 - 22^ + 25^ 



.(59) 



The ratios between the above critical times and r, 
text I t are represented in Figure 0] None of these ra- 
tios depend on x. This is explained by an interesting 
scaling property of the envelops (which is also satisfied 
separately by l^l 2 )! 



a \ipt( ax > at ) + V> s ( ax > at )\u/i 
= \iPf{x,t)+^( Xl t)\ 2 u/L 



(60) 



which holds for any fixed excitation energy and has been 
also noted in the nonrelativistic saddle contribution to 
the density [T(ij|. A consequence is that the maxima and 
minima of the envelops with respect to x and t travel at 
constant velocity, which may only depend on loq. Nev- 
ertheless, when the interference is taken into account a 
similar relation is not satisfied by Eq. i|55|) . Equation 
(|42|l implies that u S) and therefore the oscillation period 
of the density oscillations, remain constant when x and 
t are scaled, instead of increasing as the time span of the 
forerunners does. 



solution in Eq. I|44(l except in the very limit t — x. In 
the same figure we have plotted the average local instan- 
taneous frequency Lj av = —Im[{d^p/dt) /ip] versus time 
Since now tp « tpf + ip s we obtain, 



m (-2t 2 y + x 2 cos(2y)) 
r] [-t 2 - lj 2 t] 2 + {-t 2 + Ljfii] 2 ) sin(2y)] ' 



with envelops 



U)Q. 



(61) 



(62) 



The density and uj av are dephased. In particular the 
maxima of the density correspond to the minima at loq 
of the frequency, so that most of the density corresponds 
to frequencies below the cut-off. (This effect is more and 
more clear for lower injection frequencies.) This means 
that the first forerunner is essentially tunneling whereas, 
in the same time domain, the nonrelativistic wave is not. 
For the second forerunner, however, the minima of the 
density are not so close to zero, see e.g. Figure|21 so that 
there is an alternating influence of frequencies above and 
below cut-off. 



VI. EXACT SOLUTION FOR SMALL x 

The small x region may be defined by the failure of 
the inequalities in Eqs. (|48|l and 149|) . It requires in 
general an exact numerical treatment because the simple 
saddle-pole aproximation is not valid anymore. Figure 
exhibits a plot of the "density" as a function of time, 
for given values of x and loq, near rj w using Eq. 138|l . 
The exact calculation (solid line) exhibits a sharp rela- 
tivistic wavefront, reaching unity at x = t, followed by 
mild oscillations of smaller amplitude and eventually by 
the asymptotic stationary density. Figure provides also 
a comparison with the one-term approximation (dashed 
line) to Eq. I|38fl . namely, 



tp(x,t) w J (r]) 



2i— Ji(v). 
V 



(63) 



One sees that the above approximation gives an excellent 
description in the vicinity of the rclativistic wavefront. 

In Figure we have plotted the "density" versus time 
for three different values of loq. The first forerunner just 
after t = x may be characterized from the expansion of 
the two Bessel functions in Eq. I|63|l for very small rj [2(| , 



B. Instantaneous frequency 

Figure [5] shows a typical probability density for the 
first forerunner region at large x, where the two saddle 
contributions dominate clearly over the pole contribu- 
tion. We have choosen a signal frequency luq — 0.5 and 
a fixed position x — 100 (r w 710). In this scale, there 
are no differences between the exact and the aproximate 



.^0 



t + x 



(64) 



It consists of a peak that decays from 1 with slope 
— (t+x)/2 independently of luq. At variance with the non- 
relativistic equation, that immediately smoothes away 
any initial singularity, the relativistic equation propa- 
gates the sharp jump of the excitation in the input signal. 
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FIG. 6: \tp\ 2 versus time at x = 10 and uio = 0.99 for the 
exact calculation (solid line) and the first term approximation 
(dashed line). 
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FIG. 7: IV"] 2 versus time at x — 0.1 for three differents values 
of uio: uio = 0.99 (solid line), uio = 0.8 (dotted line) and 
uio = 0.1 (dotted-dashed line). The same quantity is plotted 
for the nonrelativistic solution (dashed line, uio = 0.99; long 
dashed line, uio = 0.8, and circles for uio = 0.1). 



Note in Figure that for deeper tunneling the asymp- 
totic relativistic signal becomes stronger than the non 
relativistic one: the pole contribution to the density has 
the same form in both cases |2, e _2K °, but the ko's are 
only equal at loq = 0. Otherwise, kq (nonrelativistic) 
> kq (relativistic) for < u>q < 1. 

In fact, the approximation in Eq. (|64f> is valid for 
any value of x, but its relevance is mostly appreciated at 
small x. At large x the first peak is immediately followed 
by rapid oscilations of large amplitude, so that the first 
forerunner is characterized better by the upper envelop 
in that case, except in the immediate proximity of the 
causal limit. 



0.09 1 1 1 1 1 

5 10 15 20 

t 

FIG. 8: LUav versus time at x = 0.1 for two differents values of 
the injection frequency: uio = 0.99 (solid line), luq = 0.1 (dot- 
ted line) . The same quantity is plotted for the nonrelativistic 
solution (dashed line) for luq = 0.99 and long dashed line for 
ujo = 0.1. The axis has been cut to show all curves with the 
same scale. 
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FIG. 9: time of arrival of the maximum of the lower envelop 
for the relativistic case (dots) and for the first forerunner in 
the nonrelativistic case (solid line) . uio has been choosen equal 
to 0.98. 



The average local instantaneous frequency at the 
causal limit x — t is always ujq, see Figure |Hl and Eq. 
Q64JI. which is quite different from the "high" frequen- 
cies of the nonrelativistic case, see Figure |H] (uj av (t = x; 
nonrelativistic) = 1.5 for uq = 0.99). In fact the first 
precursor as well as the rest of the wave tunnel cleanly 
at any time if the source frequency is small enough. This 
limit value of luq for total tunneling decreases when x 
increases. 

We may define as before the second forerunner ac- 
cording to some convention. Numerically it is relatively 
simple to identify the maximum of a lower envelop for 
wo > 0.8. Figure IH1 shows that, as in the nonrelativistic 
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FIG. 10: t Pmin versus using u from 0.95 to 0.999. 

case, we have a region where the signal arrives at earlier 
times at greater x, but for large x it grows linearly with 
x. The arrival time of the envelop's peak at the minimum 
of the basin, t Pmi „, depends linearly on 1/kq, see Figure 
1101 as in the nonrelativistic case. 

VII. SUMMARY AND DISCUSSION 



and position scaling as in the nonrelativistic case, so that 
qualitatively different "shallow" and "deep" tunneling 
regimes may be distinguished; (b) tunneling is more ro- 
bust relativistically, both in the precursors and asymptot- 
ically; (c) the "first" relativistic precursor, right after the 
limit imposed by causality does not have a nonrelativistic 
counterpart, and does tunnel; (d) the "second" precur- 
sor, which tends to the nonrelativistic one for excitation 
frequencies near cut-off, has an oscillating structure that 
may be characterized by its envelops. The traversal time 
t is not an exact measure of their arrival except in the 
nonrelativistic limit. 

While the emphasis here has been on the compari- 
son between the results of the Schrodinger and the rel- 
ativistic wave equation for a canonical input signal, our 
next objective is to incorporate to the theory more ele- 
ments relevant for the waveguide experiments, such as 
dissipation, other forms of input envelop pulses (e.g. 
square or Gaussian), waveguides with section constraints 
or barriers of lower dielectric constant [3^ (for classi- 
cally forbidden regions of finite width, the nonrelativis- 
tic description[ll|, |33j exhibits also the "shallow" and 
"deep" tunneling regimes mentioned above), frequency 
band- limitations, and a separate analysis of cosine or 
sine excitations, which have been combined here to relate 
directly the relativistic wave amplitudes to the nonrela- 
tivistic densities. 



In this work we have described the forerunners and the 
transition to the asymptotic regime of the solution of a 
relativistic wave equation for a unit-step-function modu- 
lated input signal with carrier frequencies below cut-off. 
This extends the investigation carried out previously for 
the Schrodinger equation 0, IToL , with the bonus that 
the relativistic equation may be physically implemented 
in waveguides. The main differences between the rela- 
tivistic and nonrelativistic cases are: (a) the relativistic 
solutions are not simply related to each other by time 
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